The number of transmission channels through a single-molecule junction 
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We calculate transmission eigenvalue distributions for Pt-benzene-Pt and Pt-butadiene-Pt junctions 
using realistic state-of-the-art many-body techniques. An effective field theory of interacting 7r-electrons 
is used to include screening and van der Waals interactions with the metal electrodes. We find that the 
number of dominant transmission channels in a molecular junction is equal to the degeneracy of the 
molecular orbital closest to the metal Fermi level. 



The transmission eigenvalues r n constitute a niesoscopic 
PIN code |0] characterizing quantum transport through any 
nanoscale device. For a single-atom contact between two 
metallic electrodes, the number of transmission channels is 
simply given by the chemical valence of the atom |2fl. Re- 
cently, highly-conductive single-molecule junctions (SMJ) 
with multiple transport channels have been formed from 
benzene molecules between Pt electrodes [3]. This raises 
the question if there exists a similarly simple criterion de- 
termining the number of transmission eigenchannels in a 
SMJ. 

Previous calculations JH 0] using effective single- 
particle models based on density functional theory appear 
to answer the above question in the negative. However, we 
find that transmission channel distributions calculated us- 
ing many-body theory do yield a simple, intuitive answer 
to this important question. 

For a two-terminal SMJ, r„ are eigenvalues of the elastic 
transmission matrix 0-01 

T(E) = 1\ (E)G(E)T 2 (E)G^ (E) , (1) 

where G is the retarded Green's function [8] of the SMJ 
and T a is the tunneling-width matrix describing the cou- 
pling of the molecule to lead a. The number of transmis- 
sion channels is equal to the rank of the matrix (1), which is 
in turn limited by the ranks of the matrices T a and G. The 
additional two-fold spin degeneracy of each resonance is 
considered implicit. The rank of T a is equal to the number 
of covalent bonds formed between the molecule and lead a. 
Thus, for example, in a SMJ where a Au electrode bonds to 
an organic molecule via a thiol group, only a single bond is 
formed, and there is only one non-negligible transmission 
channel 10,0]. In Pt-benzene-Pt junctions, however, each 
Pt electrode forms multiple bonds to the benzene molecule 

a. 

In this article, we investigate how transmission eigen- 
value distributions of SMJs depend on the number of 
lead-molecule bonds and on molecular symmetry using a 
many -body theory of transport Hfl. Specifically, we focus 
on junctions with benzene (C 6 H 6 ) and butadiene (C4H 6 ) 
bonded to two Pt leads (see Fig. Q]). Consistent with refs. 
SandB we find that the total number of nonzero transmis- 
sion eigenvalues in a SMJ is limited only by the number of 



FIG. 1. Schematic diagram of a Pt-benzene-Pt junction. The 
lead-molecule coupling is mediated predominantly by individual 
Pt atoms at the tips of each electrode. 



bonds to each electrode. However, increasing the number 
of bonds past a certain point leads to additional channels 
with very small transmission (r„) <C 1. The central find- 
ing of this article is that in SMJs with sufficient numbers of 
lead-molecule bonds the number of dominant transmission 
channels is equal to the degeneracy of the molecular or- 
bital closest to the metal Fermi level. Additional transmis- 
sion channels stemming from further off-resonant molecu- 
lar states are strongly suppressed, but may still be experi- 
mentally resolvable y|] for very strong lead-molecule hy- 
bridization. 



MANY-BODY THEORY OF TRANSPORT 

When macroscopic leads are attached to a single 
molecule, a SMJ is formed, transforming the few-body 
molecular problem into a true many-body problem. Un- 
til recently, a theory of transport in SMJs that prop- 
erly accounts for both the particle and wave character 
of the electron has been lacking, so that the Coulomb 
blockade and coherent transport regimes were considered 
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'complementary.' IToll In this article, we utilize a nonequi- 
librium many-body theory [8J] that correctly accounts for 
wave -particle duality, reproducing the key features of both 
the Coulomb blockade and coherent transport regimes. 
Previous applications 1181 [I114l3ll to SMJs utilized a semi- 
empirical Hamiltonian lfl4TI for the 7r-electrons, which ac- 
curately describes the gas-phase spectra of conjugated or- 
ganic molecules. This approach should be adequate to de- 
scribe molecules weakly coupled to metal electrodes, e.g. 
via thiol linkages. However, in junctions where the ir- 
electrons bind directly to the metal electrodes |3fl, the lead- 
molecule coupling may be so strong that the molecule itself 
is significantly altered, necessitating a more fundamental 
molecular model. To address this issue, we have developed 
an effective field theory of interacting n-electrons, in which 
the form of the molecular Hamiltonian is derived from 
symmetry principles and electromagnetic theory (multi- 
pole expansion), rather than using an ad-hoc parametriza- 
tion n X 411 - The resulting formalism constitutes a state-of- 
the-art many-body theory that provides a realistic descrip- 
tion of lead-molecule hybridization and van der Waals cou- 
pling, as well as the screening of intramolecular interac- 
tions by the metal electrodes, all of which are essential for 
a quantitative description of strongly-coupled SMJs 0]. 
The Green's function of a SMJ has the form JH] 

G(E) = [G^JE) - E T (£?) - A£ c (£)] 1 , (2) 

where G mo \ is the molecular Green's function, S T is 
the tunneling self-energy matrix, whose imaginary part is 
given by ImS T = — ^ Q T a /2, and AS C is the correc- 
tion to the Coulomb self-energy due to the broadening of 
the molecular resonances in the junction. At room tem- 
perature and for small bias voltages, AS C ~ in the co- 
tunneling regime [8] (i.e., for nonresonant transport). Fur- 
thermore, the inelastic transmission probability is negligi- 
ble compared to Eq. [TJin that limit. 

The molecular Green's function G mo \ is found by ex- 
actly diagonalizing the molecular Hamiltonian, projected 
onto a basis of relevant atomic orbitals [8J]: 

[ V {u) + V (u')]c{vy) 

^ E-E v ,+E v + iQ+ ' 



Gmol (E) 



(3) 



where E v is the eigenvalue associated with eigenstate v 
of the molecular Hamiltonian, V{y) is the probability that 
the state v is occupied, and C(u, v') is a rank-1 matrix with 
elements 

[C(u,u')] na , ma , = {v\d na \v')(v'\d\ nrJ ,\v). (4) 

Here d na annihilates an electron of spin a on the nth 
atomic orbital of the molecule. For linear response, V{y} 
is given by the grand canonical ensemble. Equations [JHH 
imply that each molecular resonance v — > v' contributes at 
most one transmission channel in Eq.[TJ 

An effective field theory of interacting 7r-electrons was 
used to model the electronic degrees of freedom most rele- 
vant for transport. Using symmetry principles and an elec- 
trostatic multipole expansion, an effective Hamiltonian for 
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FIG. 2. Energy-level diagram of a molecule with ioniza- 
tion energy IE and electron affinity EA bound to a metal sur- 
face/electrode with work function <j>. In forming the metal- 
molecule bond, the HOMO resonance of the molecule shifts and 
hybridizes. 



the 7r-electrons was derived , keeping all interaction terms 
up to the quadrupole-quadrupole interaction. The minimal 
set of parameters necessary to characterize the 7r-electron 
system consists of the 7r-orbital on-site repulsion U , the 
7r-orbital quadrupole moment Q, the nearest-neighbor tt- 
7r hybridization t, and the dielectric constant e, which 
accounts for the polarizability of the cr-electrons, which 
are not included explicitly in the calculation. The values 
C/ =8.9eV, Q = 0.67eA 2 , t = 2.64eV, and e = 1.63 
were determined by fitting to the gas-phase spectrum of 
benzene. 11411 In the molecular junction, screening of in- 
tramolecular Coulomb interactions by the nearby metal 
electrodes El was included via the image charge 
method , with no additional adjustable parameters. 



JUNCTION ENSEMBLE 

In this article, we consider junctions in which two 
macroscopic multi-channel leads each couple to several 
atomic orbitals of a single molecule. In order to calcu- 
late the distribution of transmission eigenvalues, it is first 
necessary to construct a physical ensemble of junctions. 
Both the lead-molecule coupling and the electrode geome- 
try Il7ll vary over the ensemble of junctions produced in an 
experiment B]. The lead-molecule coupling involves both 
screening lll5Lll6ll and hybridization of the molecular and 
metallic states, described by the matrix S T (see Fig. [2]). We 
assume both effects are dominated by the interaction of the 
molecule with a single Pt atom at the tip of each electrode, 
as illustrated in Fig. [TJ Since S T depends exponentially 
on the tip-molecule distance, we assume its variation is 
most important, and keep screening fixed over the ensem- 
ble of junctions. Moreover, we neglect the real part of S T 
and focus on the variation of the tunneling-width matrices 
r Q . The variation of electrode geometry leads to a varia- 
tion of the work function (j) and density of states (DOS) of 
the leads (the latter also contributes to the variation of T a .) 
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We thus assume the ensemble of junctions can be modeled 
adequately through variations of T a and (j) only 111 811 . 



Molecule-electrode Coupling 

In order to determine the tunneling-width matrices T a 
for a SMJ, we first consider the details of a single ben- 
zene molecule adsorbed on a Pt(lll) surface. This is the 
most stable Pt surface, and has been the subject of numer- 
ous investigations, where the observed binding energy for 
benzene ranges between 21 kcal/mol (0.91 eV/molecule) 
to 29 kcal/mol (1.26 eV/molecule) corresponding to the 
atop(0) and bridge(30) bonding configurations, respec- 
tively lfl9l-21l. As indicated schematically in Fig. [2 when 
a molecule binds with a metal surface, the relevant energy 
levels of the molecule shift and hybridize, forming bond- 
ing and anti-bonding states. These two effects contribute 
to the binding energy AE h = AE vdW + AE hyh , where 
AE vdW = (H mo i) — (H mo i) is the van der Waals en- 
ergy shift and AE hyh is the hybridization energy. Here 
-ff mo i is the gas-phase molecular Hamiltonian and H mo \ is 
the molecular Hamiltonian including screening from the Pt 
surface. Taking the benzene-Pt distance as 2.25A and as- 
suming the screening is dominated by the nearest Pt atom, 
we find that the HOMO-LUMO gap of benzene reduces 
from 10.05eV in the gas-phase to 7.52eV on Pt(lll) and 
A£ vdW = 0.49eV. This implies AE hyh < 0.77eV. 

Since the metallic work function <fi lies between the 
HOMO and LUMO resonances, hybridization occurs via 
the virtual exchange of an electron or hole between the 
metal and the neutral molecule. Using second-order per- 
turbation theory, we find: 



A£ hyb = 



+ E 



30 dE Tr{T{E)C{v^ N )} 
2^ E-E Qn +E v 

" dETr{T(E)C(Q N ,u')} 



2?r 



-E — E + E v 



(5) 



where fx is the chemical potential of the lead metal, ~H N 
is the iV-particle molecular Hilbert space, and N is the 
ground state of the iV-particle manifold of the neutral 
molecule. At room temperature, the Pt DOS g(E) is 
sharply peaked around the Fermi energy [22], allowing 
us to perform the energy integral in Eq. fusing T(E) ~ 
T(e F )Z5(E - e F )/g(e F \ where we take Z=+4 for a Pt 
atom and g(e F )=2.88/eV Il23ll . The hybridization energy is 
thus determined by the tunneling-width matrix evaluated at 
the Pt Fermi level, r(e F )=T. Using </> Pt(m) =5.93eV JH 
and the chemical potential of benzene /i =(IE + EA) /2=- 
4.06eV I25ll26ll (which is unaffected by screening), we find 
TV {r} < 21.6eV to fit the hybridization energy for ben- 
zene adsorbed on Pt. 




FIG. 3. Distributions of hybridization energy (top) and tunneling- 
width matrix trace (bottom) for a Pt-benzene-Pt junction. The 
99 th -percentile values for AE hyh and Tr{T Q } are 0.77eV and 
10.82eV, respectively. Here a=l,2 labels the lead-molecule con- 
tacts. 



RESULTS AND DISCUSSION 
Pt-benzene-Pt Junctions 

Let us now consider lead-molecule hybridization in a 
Pt-benzene-Pt junction. Screening from two metal elec- 
trodes further reduces the HOMO-LUMO gap to 6.46eV. 
Since the most favorable binding of benzene on the closest- 
packed Pt(lll) surface gives A£' hyb =0.77eV, we assume 
this is essentially an upper bound on hybridization in a 
SMJ, where the bonding is more random. We wish to study 
the dependence of the transmission eigenvalue distribution 
on the number of bonds formed with each electrode. For M 
covalent bonds between a macroscopic lead and a molecule 
with P atomic orbitals, L is a rank-M matrix, which can 
be represented as 



r 



M 



(6) 



where -y m are linearly-independent real row vectors of di- 
mension P, representing linear combinations of the atomic 
orbitals of the molecule [27]. Our approach is to populate 
the elements of 7 m from a uniform random distribution on 
the interval [-.A,A|. The bonding ensemble corresponds to 
a random walk of M steps in a P-dimensional space. The 
distributions of Ai^yb and Tr {T a } shown in Fig.[3]have 
long Gaussian tails, so the value A = 0.82eV was chosen 
to fix the 99 th -percentile of AE hyh at 0.77eV [51. The 
99 th -percentile of Tr {T Q } is 10.82eV which, per orbital, 
is nearly 3 x the coupling found for a Au-BDT-Au junction 

a. 

In addition to sampling a variety of bonding configu- 
rations, we assume the ensemble of junctions samples all 
possible Pt surfaces. The work function of Pt ranges from 
5.93eV to 5.12eV for the (111) and (331) surfaces, respec- 
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FIG. 4. Calculated transmission eigenvalue distributions for an 
ensemble of 24 x 10 4 Pt-benzene-Pt junctions. 3,000 bonding 
configurations x 80 Pt surfaces were included. Each lead was 
assumed to form A/=6 bonds with the molecule. There are two 
dominant transmission channels arising from the two-fold degen- 
erate HOMO resonance closest to the Pt Fermi level, with a small 
but experimentally resolvable third channel arising from further 
off-resonant tunneling. 



tively 12411 . so that 



1.88eV < /z Pt - Ho < -1.07eV, 



(7) 



and we assume a uniform distribution of fip t on this inter- 
val. 

The transmission eigenvalue distributions are shown in 
Fig- 0] for an ensemble of 24 x 10 4 Pt-benzene-Pt junc- 
tions. Each lead was assumed to form M = 6 bonds with 
the molecule. Despite the existence of six covalent bonds 
between the molecule and each lead, there are only two 
dominant transmission channels, which arise from the two- 
fold degenerate HOMO resonance closest to the Pt Fermi 
level. There is also a small but experimentally resolvable 
third channel arising from further off -resonant tunneling. 
This channel is non-negligible because of the very large 
lead-molecule coupling T in the Pt-benzene junction. Ex- 
plicit calculations with smaller values of Tr {T} yielded 
only two non-negligible transmission eigenvalues. Thus, 
for metals with a smaller DOS at the Fermi level, such as 
Cu, Ag, or Au, junctions with benzene would be expected 
to exhibit only two measurable transmission channels. 

The calculated transmission eigenvalue distribution 
shown in Fig.|4]is consistent with the experiment [ 3], which 
determined the transmission eigenvalues for three partic- 
ular junctions: {r„}={0.68, 0.40}, {0.36, 0.25, 0.10}, 
{0.20}, where a third channel was observable only once. 
The conductance histogram for the same ensemble of junc- 
tions is shown in Fig. [5] The peak conductance value is 




0.25 0.5 0.75 1 1.25 1.5 1.75 2 

G (2e 2 /h) 

FIG. 5. Calculated conductance histogram for the same en- 
semble of Pt-benzene-Pt junctions discussed in Fig. |4j where 

G=(2e 2 /h)j: n Tn. 



some 20% less than that reported experimentally [3]. This 
discrepancy might be attributable, in part, to the inclusion 
of a small fraction of Pt-Pt junctions in the experimental 
histogram. 

In addition to the ensemble of junctions shown in Fig. SI 
we also investigated ensembles of junctions with M a = 
1,...,6, including the case Mi 7^ M 2 . Consistent 
with the discussion in Refs. El, we find that the total 
number of nonzero transmission eigenvalues is M min = 
min{Mi, M 2 }. However, whenever M min > 2 there are 
always two dominant transmission channels, and the total 
transmission probability does not increase appreciably be- 
yond M min = 2. 

The above analysis demonstrates that the two dominant 
transmission channels evolve from the two-fold degener- 
ate HOMO resonance in Eq.[3]as the lead-molecule cou- 
pling S T is turned on. For finite S T , the poles of G mo i 
are mixed by Dyson's equation (Eq.|2|, making it problem- 
atic to decompose the transmission eigenvalues into sepa- 
rate contributions from each molecular resonance fH]. Al- 
ternatively, the projections of the transmission eigenvec- 
tors onto the molecular resonances can be computed |4J]. 
Because an "extended molecule" is often used in density- 
functional calculations to account for charge transfer be- 
tween the molecule and electrodes, it is difficult if not 
impossible to interpret these contributions in terms of the 
resonances of the molecule itself [4]. Since charging ef- 
fects in SMJs are well-described in our many-body theory 
fH], there is no need to utilize an "extended molecule," so 
the projections of the transmission eigenvectors onto the 
molecular resonances can be determined unambiguously 
(see Methods). We find that the mean-square projections 
of the first and second transmission channels in Fig.|4]onto 
the benzene HOMO resonance are 87% and 71%, respec- 
tively, confirming the conclusion that these eigenchannels 
correspond to tunneling primarily through the HOMO res- 
onance. 



Pt-butadiene-Pt Junctions 

To test our hypothesis that the number of dominant 
transmission channels is limited by the degeneracy of the 
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FIG. 6. Transmission eigenvalue distributions for an ensemble 
of 24 x 10 4 Pt-butadiene-Pt junctions. 3,000 bonding configu- 
rations x 80 Pt surfaces were included. Each lead was assumed 
to form M = 4 bonds with the molecule. The linear butadiene 
(C4H6) molecule lacks the orbital degeneracy of the benzene ring 
and consequently exhibits only 1 dominant channel. A second 
channel due to further off-resonant tunneling may be experimen- 
tally resolvable. The total conductance distribution peaks around 
~0.4G . 



resonance, which lies closest to the Pt Fermi level, result- 
ing in two dominant transmission channels for benzene 
(2-fold degenerate HOMO) and a single dominant trans- 
mission channel for butadiene (non-degenerate HOMO). 
Our results for the transmission channel distributions of Pt- 
benzene-Pt junctions (see Fig. |4j are in quantitative agree- 
ment with experiment yfl. 

Despite the larger number of states available for tun- 
neling transport in SMJs, we predict that the number of 
transmission channels is typically more limited than in 
single-atom contacts because molecules are less symmet- 
rical than atoms. Nonetheless, certain highly-symmetric 
molecules exist that should permit several dominant trans- 
mission channels. For example, the C 60 molecule pos- 
sesses icosahedral symmetry, and has a 5-fold degener- 
ate HOMO resonance and 3-fold degenerate LUMO res- 
onance Bill . C 6 o-based SMJs with gold electrodes have 
been fabricated, and shown to exhibit fascinating elec- 
tromechanical 02I1 and spin-dependent [33] transport prop- 
erties. For Pt-C 60 -Pt junctions, we predictive dominant 
transmission channels stemming from the C 60 HOMO res- 
onance, which lies closest to the Pt Fermi level [34ll35ll. 



most relevant molecular orbital, we have investigated trans- 
mission eigenvalue distributions for Pt-butadiene-Pt junc- 
tions. Butadiene is a linear conjugated polymer, lacking the 
discrete (six-fold) rotational symmetry of benzene. Since 
the molecular orbitals of butadiene are nondegenerate, we 
anticipate a single dominant transmission channel. Us- 
ing the same per-orbital hybridization as for benzene gives 
TV{r a } < 7.21eV. The image charge method gives 
A£ vdW = 0.52eV and we find AE hyh < 0.59eV, so 
that AE h < 1.12eV = 107.7kJ/mol, in line with exist- 
ing values reported in the literature l30h . The range of Pt 
work functions for all possible Pt surfaces gives a chemical 
potential range of 

- 1.70eV < fipt - fa < -0.89eV. (8) 

Despite forming four bonds with each electrode, it is ev- 
ident from Fig. [6] that the Pt-butadiene-Pt junction has a 
single dominant transmission channel. The mean-square 
projection of this transmission channel onto the non- 
degenerate butadiene HOMO resonance is 80% (see Meth- 
ods). 

CONCLUSIONS 

In conclusion, we find that the number of dominant 
transmission channels in a SMJ is equal to the degener- 
acy of the molecular orbital closest to the metal Fermi 
level. Transmission eigenvalue distributions were calcu- 
lated for Pt-benzene-Pt and Pt-butadiene-Pt junctions us- 
ing realistic state-of-the-art many-body techniques. In both 
cases, transmission occurs primarily through the HOMO 



METHODS 

Decomposing transmission channels in terms of molecular 
resonances 

In order to understand how transport in a SMJ is deter- 
mined by the chemical properties of the molecule, it would 
be desirable to express the transmission eigenchannels in 
terms of molecular resonances 0, 0]. The transmission 
eigenvalues r„ and eigenvectors jn) are solutions of the 
equation 

T\n)= Tn \n), (9) 
where T is the transmission matrix given by Eq.[TJ By def- 
inition, \n) is a linear combination of the atomic orbitals of 
the molecule. In an effective single-particle model [4J], |n) 
may also be expressed as a linear combination of molecular 
orbitals \<pj) 

\n) = J2 a M- ( 10 ) 

3 

Here \ct£\ 2 can be identified as the contribution of the 
jth molecular orbital to the nth transmission channel |4J], 
which can be conveniently expressed in terms of the pro- 
jection operator P 3 = \4>j){4>j \ as \a£\ 2 = (n\Pj\n). 

In the many-body problem, there is no orthonormal set 
of "molecular orbitals"; rather each molecular resonance 
of energy E v > — E v corresponds to a transition v — > v' be- 
tween an iV-body and an (N + l)-body molecular eigen- 
state [see Eqs. [3] and [4). The projection operator onto a 
molecular resonance is 

* C(v, 1/) 

v ^ v ~ Tr{G>,z/)}' K ' 
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where C(u, v') is given by Eq.|4l The absolute square pro- 
jection of the nth transmission eigenvector onto the reso- 
nance v — > v' is given by 



a,, 



\ z = {n\P u ^ v ,\n). (12) 

A necessary condition to identify an eigenchannel \n) with 
transmission through a particular molecular resonance is 

I u^v'\2 ^ l 

The above procedure is in principle straightforward to 
implement in an effective single-particle model based on 
density functional theory (DFT)||4j]. However, in practice, 
an "extended molecule" must be used in DFT calculations 
to account for charge transfer between the molecule and 
electrodes. This is because current implementations of 
DFT fail to account for the particle aspect of the electron 
ESS, i.e., the strong tendency for the electric charge 
on the molecule within the junction to be quantized in inte- 
ger multiples of the electron charge e. Analyzing transport 
in terms of extended molecular orbitals has unfortunately 
proven problematic. For example, the resonances of the 
extended molecule in Ref. 4 apparently accounted for less 
than 9% of the current through the junction. 

Since charging effects in SMJs are well-described in our 
many -body theory JD], there is no need to utilize an "ex- 
tended molecule," so the projections of the transmission 
eigenvectors onto the molecular resonances can be deter- 
mined directly from Eq.[T2l 



Benzene resonances 

The neutral ground state of benzene is nondegenerate, 
while the HOMO and LUMO resonances are both doubly 
degenerate due to the (six-fold) rotational symmetry of the 
molecule. To be consistent with the discussion of Ref. |j, 
the additional two-fold spin degeneracy of each resonance 
is considered implicit. 

We define the following projection operators: 



Pi 



HOMO 



E 

^eo 5 



P, 



v— >0 6 ) 



P 



LUMO 



= TP t 



Oo- 



Pl 



HOMO 



P 



LUMO; 



(13) 



(14) 



(15) 



where N is the ground state with N 7r-electrons and 
1 is the six-dimensional unit matrix in the space of ir- 

orbitals. Promo an d Plumo are projection operators onto 
the two-dimensional subspaces spanned by the HOMO 
and LUMO resonances, respectively. Eqs. 1X314131 define 
a complete, orthogonal set of projection operators in six- 
dimensions with the properties 



E 



Pi, Pi Pi 0« 



(16) 



(17) 



In particular, Eq. [171 implies that the HOMO and LUMO 
subspaces of benzene are orthogonal. The absolute square 
of the projection of the n th transmission eigenvector onto 
the subspace spanned by P, is 



>\* = (n\Pj\n). 



(18) 



These coefficients satisfy the condition 



(19) 



where the sum runs over j=HOMO, LUMO, _L. 

Fig.|7]shows the mean-square projections (|o^ | 2 ) of the 
transmission eigenvectors onto (a) the benzene HOMO 
resonance; (b) the benzene LUMO resonance; and (c) the 
two-dimensional subspace orthogonal to both the HOMO 
and LUMO resonances, as a function of electrode chemi- 
cal potential for the same ensemble of Pt-benzene-Pt junc- 
tions discussed in Fig. [4] We find that the mean-square pro- 
jections of the first and second transmission channels onto 
the benzene HOMO resonance are 87% and 71%, respec- 
tively, confirming the conclusion that these eigenchannels 
correspond to tunneling primarily through the HOMO res- 
onance. 

Midway between the HOMO and LUMO resonances 
at u = fj, , the first two transmission channels have 
mean-square projections of ~ 0.5 onto both the HOMO 
and LUMO resonances, consistent with the expectation 
that the HOMO and LUMO resonances should contribute 
equally to transmission. The remaining channels do not 
have negligible overlap with the HOMO resonance, but 
instead cluster around (| a^ OMO 1 2 ) ~ 0.3. The transmis- 
sion channels with (r n ) <C 1 correspond to contributions 
from several far off -resonant poles, each of which has some 
overlap with the HOMO resonance due to the overcom- 
pleteness of the projectors P u ^, v > [39]. They are essentially 
random unit vectors in the six-dimensional space of ben- 
zene 7r-orbitals, whose mean-square overlap with the two- 
dimensional HOMO subspace should be (|a^ OMO | 2 ) = 
2/6 = 1/3. 

The first two transmission channels only have an appre- 
ciable overlap with the subspace orthogonal to the HOMO 
and LUMO resonances in the vicinity of a pronounced dip 
in the transmission spectrum at u — u ~ ±3.5eV. 
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FIG. 7. The mean-square projections |an| 2 = (n\Pj\n) of the 
transmission eigenvectors onto the two-dimensional subspaces 
(a) spanned by the benzene HOMO resonance; (b) spanned by 
the benzene LUMO resonance; and (c) orthogonal to the HOMO 
and LUMO subspaces, for the same ensemble of Pt-benzene-Pt 
junctions discussed in Fig. [4] For the range of possible chemical 
potentials of Pt electrodes, -1.88eV< fip t — fiQ < -1.07eV indi- 
cated by the grey boxes in each subfigure, the first two channels 
have very strong overlap with the HOMO resonance: 87% and 
71%, respectively. 



Butadiene resonances 

The neutral ground state of butadiene is nondegenerate, 
and the HOMO and LUMO resonances have no orbital 
degeneracy. The projection operators onto the HOMO and 
LUMO resonances of butadiene are 



Phomo 
Plumo 



Pq 3 



Pr 



4 -s-0 5 ) 



(20) 
(21) 



respectively. Fig. [8] shows the mean-square projections of 
the transmission eigenvectors onto the nondegenerate bu- 
tadiene HOMO and LUMO resonances as a function of 
electrode chemical potential for the same ensemble of Pt- 
butadiene-Pt junctions discussed in Fig. [6] In this case, the 
single dominant channel has a strong overlap with the non- 
degenerate HOMO resonance: (|af OMO | 2 ) = 0.80 when 
averaged over the range -1.70eV< /Up t — Ho < -0.89eV, 
while the far off -resonant channels with (r n ) <C 1 have 
(|a^ OMO | 2 ) ~ 1/4, as expected based on the arguments 
given above. 
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